Dark— Bright Solitons in Inhomogeneous Bose-Einstein Condensates 
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We investigate dark-bright vector solitary wave solutions to the coupled non-linear Schrodinger 
equations which describe an inhomogeneous two-species Bose-Einstein condensate. While these 
structures are well known in non-linear fiber optics, we show that spatial inhomogeneity strongly 
affects their motion, stability, and interaction, and that current technology suffices for their creation 
and control in ultracold trapped gases. The effects of controllably different interparticle scattering 
lengths, and stability against three-dimensional deformations, are also examined. 



Among the many surprising features of non-linear 
equations, the emergence of solitons is one of the most 
prominent. For the non-linear Schrodinger equation 
(NLSE), which governs both non-linear optical modes 
in fibers and dilute Bose-Einstein condensates (BECs), 
two different kinds of scalar solitons, bright and dark, 
are known Q. In optics, bright and dark solitons arise 
in media with anomalous dispersion and normal disper- 
sion, respectively, and for BECs the s-wave scattering in- 
teraction is the determining factor (attractive for bright 
solitons, repulsive for dark). Whereas in gaseous Bose- 
condensates dark solitons only, and only recently, have 
been observed fi [| , optical solitons are well investigated 
and on the verge of industrial application [Q . In addition 
to the bright and dark scalar solitons there are also vari- 
ous multicomponent (vector) solitons known, which arise 
as solutions to systems of coupled NLSEs. An elegant ex- 
ample is the so-called dark-bright soliton, where a bright 
optical solitary wave exists in a system with normal dis- 
persion because it is trapped within a copropagating dark 
soliton H ||, 0] . In this Letter we investigate the behav- 
ior of dark-bright solitons (and solitary waves) in repul- 
sively interacting two-component Bose-Einstein conden- 
sates. We examine the effects of spatial inhomogeneity, 
three dimensional geometry, and dissipation, which are 
all important features of BEC experiments. 



FIG. 1: (a) A dark-bright soliton solution of eqs. (^), with 
a — 0. The rescaled densities of the bright and dark compo- 
nents, \4>b\ 2 and IV'dI 2 , are shown with a broken and a full 
line, respectively. (6) The size of a motionless dark-bright 
soliton, in units of the healing length /i -1//2 , as a function of 



ural units, the general equations read 

ii>B = -^b + [Vb + \iPb\ 2 +Gd\tPd\ 2 -V- A]Vb , 

ii> D = -\i>D + Wd + \iPd\ 2 + g b \M 2 - (A^d , (i) 

where the chemical potentials /id — M and [Lb = M + 
A have been introduced in the standard way. For the 
present we set the coupling strength ratios Gb,d = 1 
and assume no external potentials, Vd,b = 0- The dark- 
bright soliton solution to eqs. (Q) is then given by 



In the context of cold atomic gases the two vector 
components evolving under the Gross-Pitaevskii NLSE 
are the macroscopic wave functions of Bose-condensed 
atoms in two different internal states, which we will de- 
note as \D) and \B). The non-linear interactions are 
due to elastic s-wave scattering among the atoms, and 
are effectively repulsive (positive scattering length) for 
both systems ( 23 Na and 87 Rb) in which multicomponent 
condensates have been realized. In both cases the three 
ID interaction strengths (<?dd, <?bb, 9db) may easily be 
made equal to within a few percent; and quasi-one dimen- 
sional traps that are longitudinally very flat are under ac- 
tive experimental development Q . So we will begin with 
the two-component NLSE with all interaction strengths 
equal, in one dimension, with no potential, and then add 
more realistic features in succession. In the standard nat- 
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■0D = iy/Jl sin a + y/Jl cos a ta,nh(^K[x — q(t)] J , (2) 

where Nb = Jdx\ipB\ 2 is the rescaled number of par- 
ticles in state \B), the soliton inverse length is K = 
^/^cos 2 a + (Nb/^) 2 — As/4, the bright component fre- 
quency shift is VIb = k 2 (1 — tan 2 a)/ 2— A, and the soliton 
position is q(t) = q(0) + ifttana. The 'binding energy' 
of the bright component in the well formed by the yb D 
mean field, in the co-moving frame, is clearly k 2 /2; the 
bright component phase shift ip is only of significance if 
there are two or more solitons. Readers familiar with the 
scalar solitons of the one-component NLSE will recognize 
ipD as a dark (or 'grey') soliton of velocity- angle a, and 
ipB as a bright soliton, which can only be found in single- 
component condensates if they have negative scattering 
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length (and so are prone to collapse) (see Fig. 0(a)). The 
Thomas-Fermi-like expansion with N B of the trapped 
bright component makes the soliton size k _1 longer than 
for a single-component dark soliton at the same /x (see 
Fig- 0(b)). 

The integrable system of eqs. (0) with Gd.b = 1 and 
Vd.b = 0, also known as the Manakov equation, con- 
serves the free energy 
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Since G decreases with increasing soliton velocity, the 
soliton is formally unstable (to acceleration!). But 
one implication of integrability is that perturbations of 
eqs. (||) due to interactions with other waves (solitary 
or ordinary) will not cause dissipation. If an inhomo- 
geneous potential is added, however, by allowing non- 
zero Vd,b(x) in (0), then the system is no longer in- 
tegrable, and the soliton can interact non-trivially with 
the surrounding condensate. Nevertheless, if V varies 
slowly on the soliton scale n, then excitations of the 
background that have high enough temporal frequency 
to accept energy resonantly from the soliton must also 
have short enough spatial wavelength to 'see' Vd,b{x) 
as approximately constant, and hence tend to decou- 
ple from the soliton as they do in the truly integrable 
case. The result is that G(q, q) (given by replacing 
H -» n - V D (q),A — > A — V B (q) + V D (qj in eq. (g) 
and inverting q — n(q, cos a) tan a) is approximately con- 
served, and this determines the motion of the soliton in 
the potential. (A more sophisticated multiple time scale 
boundary layer analysis supports this simple argument.) 

The motion this implies simplifies in the limit of veloc- 
ities much smaller than the speed of sound, because 
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dropping a term which is constant if we assume that none 
of the conserved state \B) atoms escape from the soliton. 
This implies the low-velocity equation of motion 



Vp(g) N B [V D (q)-2V B (q)} 
2 Sy/n + N- B l\Q-V D {q) 



(5) 



which, together with its numerical confirmation shown in 
Fig. 0, is the primary result of this paper. In the limit 
N B — ► we recover the equation of motion of the dark 
soliton || , and as N B increases we find that the soliton is 
more and more insulated from the effect of Vd , and more 
sensitive to V B — Vd- In the limit N B 3> *//J, where 
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FIG. 2: Period of oscillation for a dark-bright soliton 
in a harmonic trap calculated numerically (diamonds) and 
from eq. (g) (solid line). The trapping potentials for 
the different components are related as Vb = "/Vd, with 
7 = —1, —0.5, 0.5, 0.75, 1, 1.25, 1.5 respectively for the graphs 
starting from below. The normalization of the dark compo- 
nent is J dx\ipD\ 2 = 1000 in all cases, so that fi ranges from 
66.2 to 67.1. The divergence of the curves for 7 = 1.25 and 
7 = 1.5 shows the breakdown of the oscillations as explained 
in the text. 



the soliton is expanded by the large bright component to 
many healing lengths in size, we have 
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so that a small differential force on the bright component 
will predominate. Our assumption that the whole soliton 
is small compared to the trap scale, however, means that 
the dark component retains its dramatic effect of giving 
the soliton an effectively negative mass: the soliton accel- 
erates in the opposite direction to a force exerted through 
Vb- If Vb and Vd are kept equal, on the other hand, 
a highly expanded dark-bright soliton with N B ^JJL 
moves in the potential as if it had a very large positive 
mass (because as one can see from eq. (0) the soliton's 
potential energy is also ~ — Vd)- Numerical integration 
of the coupled NLSEs shows excellent agreement with 
eq. (5) (see Fig. ||). Note that for harmonic V B > Vd 
eq. (5) implies that there is a critical Nb above which 
the soliton will escape from the trap instead of oscillat- 
ing. While the precise transition point between very slow 
oscillation and very slow escape is difficult to check nu- 
merically, our numerical results confirm that escape does 
occur in this case at larger N B . 

A trapping potential also modifies the interactions be- 
tween solitons. While generally solitons that are more 
than a few soliton lengths apart are essentially unaffected 
by each other, at closer ranges they can distort each other 
significantly. Although the respective bright component 
numbers oc N B of two solitons are simply conserved dur- 
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FIG. 3: Symmetric collision of two dark-bright solitons in a 
harmonic trap (7 = 1). Degree of brightness indicates \i/jb\ 2 
as a function of x and t (both in trap units), for repulsive 
(A<f> = 0) and attractive (A0 = n) interaction. 





FIG. 4: Collision of two dark-bright solitons in a harmonic 
trap (7 = 1). Shown is |i/>s| 2 as a function of x and t, both 
in trap units. Initially one soliton is at rest (i.e. q = 0); after 
the collision both solitons are oscillating. 



ing such interactions, the relative phase of the two bright 
components strongly affects the details of the interac- 
tion |h], : bright-dark solitons repel each other when 
the phase difference between the bright components is 
A<p = <j)\ — <j>2 = 0, and attract each other when this 
difference is Acj) = ir. This short range behaviour, which 
is opposite to that of scalar bright solitons, occurs inde- 
pendently of the confining potential, but if the effect of 
a potential is to keep two solitons close together, then 
their phase-dependent interaction can significantly affect 
their oscillations: see Fig. (|3|). Finally, even if the trap 
does not confine both solitons within their interaction 
range, the inhomogeneous potential will still qualitatively 
modify the effects of soliton collisions, fn the integrable 
case with no potential, soliton collisions are trivial, in 
the sense that the asymptotic states that emerge after 
a collision are the same as they were before it. In gen- 
eral, however, there is a net effect of even an integrable 
collision: a 'jump-like' spatial translation of the solitons, 
relative to where each would have been if it had not en- 
countered the other. Fig. || shows that in a trap such a 
translation can imply a transfer of energy between soli- 
tons as a result of a collision. 

In addition to the inhomogeneous potential, integrabil- 
ity may also be destroyed under experimental conditions 
by the fact that the three interaction strengths gDD,9BB 
and goB governing the non-linearity will generally dif- 
fer. In the quasi-ID limit, the interaction strengths are 
given by gij = / (A, + Aj), where is the 3D s-wave 
scattering length and Aj is the cross-sectional area of 



the trap confining species j = D, B. Since in our stan- 
dard natural units we have absorbed in i/jd,b the self- 
interaction strengths goD and gsB, in this more general 
case we must allow the co-efficients Gb = gDB / gBB and 
Gd = gDB / gDD to differ from unity in the NLSE system 
(|l|) . It is straightforward to show that 
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so that varying the relative tightness of radial confine- 
ment for the two species yields one free control param- 
eter, which can allow significant retuning of Gd,b even 
without the measure of modifying the scattering lengths 
themselves. 

If Gd,b 7^ 1; solutions that are distorted versions of the 
dark-bright soliton certainly exist, although they may 
only be given in closed form for special cases. (For in- 
stance, in the limit of small Nb one may discard the 
\iPb\ 2 terms in the NLSE, and find dark soliton solutions 
for tpD, with ipB oc sech"[/«(x — q)] for v(y + 1) = 2Gd- 
For Gd > 1, there are also one or more excited bound 
states oitpB-) Such solutions are often referred to as soli- 
tary waves, rather than solitons, to indicate that they 
may interact nontrivially with other solitary or ordinary 
waves. This means, for example, that a collision between 
dark-bright solitary waves may effect a net transfer of 
bright component from one solitary wave to the other; 
see Fig. (||). It also implies that unlike true solitons, 
which are transparent to all quasiparticle modes, dark- 
bright solitary waves will suffer from dissipation due to 
collisions with thermal particles and phonons, even when 
the one-dimensional approximation is excellent. It is a 
problem beyond the scope of this Letter to compute scat- 
tering rates with Gd.b significantly different from unity. 
For Gd.b close to unity, however, simple estimates for 
the anti-damping rate due this effect show it to be negli- 
gible, at attainably low temperatures, because the cross 
sections for dissipative collisions are proportional to the 
squares of the scattering length differences. In current 
experiments, however, soliton lifetimes are limited not 
by one-dimensional dissipation, but by the breakdown of 
the one-dimensional approximation. In more than one di- 
mension, single-component dark solitons suffer from the 
well-known 'snake mode' dynamical instability, in which 
bending of the plane of the density minimum grows expo- 
nentially even without dissipation. This is the presumed 
cause of decay of solitons that have been produced in 
BEC experiments to date. As shown by Muryshev et 
al. |E|, a single-component dark soliton is only unsta- 
ble to snake modes of wavelength greater than the soliton 
size, so that radial confinement to within a healing length 
should stabilize dark solitons. By an extension of their 
method of analysis, one can show that the dark-bright 
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FIG. 5: Collisions between an initially empty (Nb = 0) dark 
soliton and a dark-bright solitary wave initially at rest in a 
harmonic trap. Each row is a separate evolution, with the left 
and right plots in each row showing \ipo\ 2 and \iPb\ 2 respec- 
tively. Horizontal axes are time, vertical axes space, both in 
trap units. In all cases J dx\ipo\ 2 = 400 and J dx\tpB\ 2 = 4, for 
\i = 36 and Thomas-Fermi radius — - 8.6. For each row, from 
top to bottom, we have (Gd, Gb) as follows: (a): (1.03,1.03); 
this case is not distinguishable from (1,1). (b): (1.5,1.5); 
transfer of bright component occurs, (c): (0.5,0.5); scatter- 
ing as well as transfer of bright component is seen in the last 
collision, (d): (0.5,1.5); viewing a true 'movie' of \tPd\ 2 re- 
veals that, during the collisions, the background cloud is much 
more significantly excited in this case than in the others; and 
this is the reason for the noticeably different soliton motion 
in this case. The case (1.5,0.5), which is not shown here, is 
just noticeably different from (1.5,1.5). 



soliton is also stable against snake modes of wavelength 
less than its size k . Since for Nb S> y/fl this can 
be much larger than the healing length, stabilization of 
dark-bright solitons against snaking should easily be pos- 
sible in current traps. (The proof is simple in the case of a 
motionless dark-bright soliton in three-dimensional bulk. 
We take eqs. (0) for a = 0, and expand the free energy G 
to second order in Stpj = 3>_y (x) cos k(y cos 9 + z sin 9 — (3) , 
where Rj and Sj are real, and j = B, D. The result for 
this perturbation of wave number k is 
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dxTk(x) / dydz cos 2 k(y cos 9 + z sine) — f3) 



shows that the spectrum of this particular Hk is positive 
for |fc| > k.) 

Even more conveniently, we note that a robust method 
for the controlled creation of dark-bright solitons has al- 
ready been presented and analysed in detail (without be- 
ing explicitly recognized as such) Jl3| . Dum et al. have 
shown that dark solitons may be created in one conden- 
sate component by adiabatic transfer of population from 
a condensate in another internal state, and have already 
noted that in the late stages of this procedure the second 
component appears as a stretched dark soliton around 
the remaining population of the first, whose wave func- 
tion approaches a hyperbolic secant. Stopping short of 
complete adiabatic transfer will in fact produce a dark- 
bright soliton of arbitrary Nb- Engineering by masked 
Rabi transfer with phase imprinting may also be possi- 
ble, and the smoother total density profile and larger size 
of a highly-expanded dark-bright soliton should simplify 
the creation of slow and stable solitons by this method. 

Our conclusion therefore is that dark-bright solitons 
move in trapped condensates much as dark solitons 
(though more slowly, if Vp = Vb), and have strong ad- 
vantages in stability and controllability. Since \iI>b\ 2 can 
be imaged separately and non-destructively, they also of- 
fer a realization of Reinhardt's and Clark's proposal to 
track solitons by trapping distinguishable atoms inside 
them[0. And in addition to advancing soliton studies 
into the inhomogeneous regime of BECs, production of 
dark-bright solitons in BECs would be the development 
of atom optical tweezers with potentially sub-micron pre- 
cision: the trapping and manipulation of ultracold atoms 
by ultracold atoms. 
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Hence, for the k sector to have positive definite free en- 
ergy, it is obviously a sufficient condition that Hk have a 
positive spectrum. Considering the possibility of purely 
imaginary <E>£> shows that this is also a necessary con- 
dition. Considering the eigenfunction sech kx (which is 
clearly the ground state because it is real and nodeless) 
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